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Abstract 

Subjet distributions were measured in neutral current deep inelastic ep scattering 
with the ZEUS detector at HERA using an integrated luminosity of 81.7 pb~ x . 
Jets were identified using the kj- cluster algorithm in the laboratory frame. Sub- 
jets were defined as jet-like substructures identified by a reapplication of the 
cluster algorithm at a smaller value of the resolution parameter y cut . Measure- 
ments of subjet distributions for jets with exactly two subjets for y cut = 0.05 are 
presented as functions of observables sensitive to the pattern of parton radiation 
and to the colour coherence between the initial and final states. Perturbative 
QCD predictions give an adequate description of the data. 
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1 Introduction 



Jet production in ep collisions provides a wide testing ground of perturbative QCD 
(pQCD). Measurements of differential cross sections for jet production [1-3] have allowed 
detailed studies of parton dynamics, tests of the proton and photon parton distribution 
functions (PDFs) as well as precise determinations of the strong coupling constant, a s . 

Gluon emission from primary quarks was investigated [4, 5] by means of the internal 
structure of jets; these type of studies gave insight into the transition between a parton 
produced in a hard process and the experimentally observable jet of hadrons. The pattern 
of parton radiation within a jet is dictated in QCD by the splitting functions. These 
functions, P a b{z, //) with a,b = q or g, are interpreted as the probability that a parton of 
type b, having radiated a parton of type a, is left with a fraction z of the longitudinal 
momentum of the parent parton and a transverse momentum squared smaller than /i 2 , 
where \i is the typical hard scale of the process. The splitting functions are calculable as 
power series in a s . Thus, the characteristics of jet substructure provide direct access to 
the QCD splitting functions and their dependence on the scale. 

The understanding of jet substructure is also important in the context of jet identification 
in boosted systems, like hadronic top decays [6] or bb final states at LHC [7]. The first 
example calls for a direct application of jet substructure, the second requires knowledge 
about jet substructure to distinguish between single- and double-quark induced jets. This 
paper presents a study of jet substructure in a more controlled hadronic-type environment 
than that provided by hadron-hadron colliders. 

Jet production in neutral current (NC) deep inelastic scattering (DIS) was previously used 
to study the mean subjet multiplicity [4] and the mean integrated jet shape [5] with values 
of a s (Mz) extracted from those measurements. In the present study, the pattern of QCD 
radiation is investigated by means of the subjet topology, providing a more stringent test 
of the pQCD calculations. 

In this paper, measurements of normalised differential subjet cross sections for those jets 
which contain two subjets at a given resolution scale are presented. The measurements 
were done as functions of the ratio between the subjet transverse energy and that of the 
jet, E^^/E^ 1 , the difference between the subjet pseudorapidity 1 (azimuth) and that of 
the jet, 7f b j — rf ct (|</> sb j — <^ ct |), and a sh \ the angle, as viewed from the jet centre, between 
the subjet with higher transverse energy and the proton beam line in the pseudorapidity- 



The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity is 
defined as t] = — ln(tan |), where the polar angle 9 is taken with respect to the proton beam direction. 



1 



azimuth plane (see Fig. 1). The predictions of pQCD at next-to- leading order (NLO) 
were compared to the data. 

2 Jets and subjets 

The analysis of subjets presented in this paper was performed using the laboratory frame. 
In this frame, the calculations of the subjet distributions can be performed up to (D(a 2 ), 
i.e. NLO, with jets consisting of up to three partons. The analysis used events with high 
virtuality of the exchanged boson, Q 2 ; at low values of Q 2 , the sample of events with at 
least one jet of high Erj? (E^ 3> y/Q 2 ) is dominated by dijet events. In that case, the 
calculations include jets consisting of up to only two partons and, therefore, correspond 
to lowest-order predictions of jet substructure. 

The kx cluster algorithm [8] was used in the longitudinally invariant inclusive mode [9] to 
define jets in the hadronic final state. Subjets [10] were resolved within a jet by consid- 
ering all particles associated with the jet and repeating the application of the kx cluster 
algorithm until, for every pair of particles i and j the quantity dij = mm(ET,i, Etj) 2 • 
((Vi ~~ Vj) 2 + (0i — 0j) 2 )> where Er,i, f]i an d 4>i are the transverse energy, pseudorapidity 
and azimuth of particle i, respectively, was greater than d cnt = y cut ■ (E^) 2 . All remaining 
clusters were called subjets. 

The subjet multiplicity depends upon the value chosen for the resolution parameter y cnt . 
Subjet distributions were studied for those jets with exactly two subjets at a value of 
the resolution parameter of y cut = 0.05. This value of y cut was chosen as a compromise 
between resolution, size of the hadronisation correction factors and statistics. The effect 
of the parton-to-hadron corrections on the shape of the subjet distributions becomes 
increasingly larger as y cnt decreases. On the other hand, the number of jets with exactly 
two subjets decreases rapidly as y cnt increases. 

Subjet distributions were studied as functions of E^ /E^, r) sb i —rf et , \<j) sh i — <ft ct \ and a sb K 
One of the goals of this study was to investigate the extent to which pQCD calculations are 
able to reproduce the observed distributions. In addition, the dependence of the splitting 
functions P a b(z, /-0 on z can be investigated using the E^/E^ distribution. The splitting 
functions at leading order (LO) do not depend on \x but acquire a weak dependence due to 
higher-order corrections. Such a dependence can be investigated by measuring the subjet 
distributions in different regions of E^ or Q 2 . 

The substructure of jets consisting of a quark-gluon pair (the quark-induced process eq — > 
eqg) or a quark-antiquark pair (the gluon-induced process eg — > eqq) are predicted to 
be different (see Section 8.1). Furthermore, the relative contributions of quark- and 
gluon-induced processes vary with Bjorken x and Q 2 . The predicted difference mentioned 
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above is amenable to experimental investigation by comparing the shape of the subjet 
distributions in different regions of x and Q 2 . 

Colour coherence leads to a suppression of soft-gluon radiation in certain regions of phase 
space. The effects of colour coherence between the initial and final states have been 
studied in hadron-hadron collisions [11]. These effects are also expected to appear in 
lepton-hadron collisions. For the process eq — > eqg, colour coherence implies a tendency 
of the subjet with lower (higher) transverse energy, i^iow (^rjiigh): ^° nave V sh * —if ct > 
(?7 sb j — rf et < 0). The variable a sh i, defined in close analogy to the variables used to study 
colour coherence in hadron-hadron collisions [11], reflects directly whether the subjet 
with the lower transverse energy has a tendency to be emitted towards the proton beam 
direction. 



3 Experimental set-up 

A detailed description of the ZEUS detector can be found elsewhere [12, 13]. A brief 
outline of the components most relevant for this analysis is given below. 

Charged particles were tracked in the central tracking detector (CTD) [14], which operated 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consisted of 72 cylindrical drift-chamber layers, organised in nine superlayers covering the 
polar- angle region 15° < 6 < 164°. The transverse-momentum resolution for full-length 
tracks can be parameterised as <t(pt)/pt = 0.0058pr © 0.0065 © 0.0014/p T , with p T in 
GeV. The tracking system was used to measure the interaction vertex with a typical 
resolution along (transverse to) the beam direction of 0.4 (0.1) cm and to cross-check the 
energy scale of the calorimeter. 

The high-resolution uranium-scintillator calorimeter (CAL) [15] covered 99.7% of the 
total solid angle and consisted of three parts: the forward (FCAL), the barrel (BCAL) 
and the rear (RCAL) calorimeters. Each part was subdivided transversely into towers and 
longitudinally into one electromagnetic section and either one (in RCAL) or two (in BCAL 
and FCAL) hadronic sections. The smallest subdivision of the calorimeter was called a 
cell. Under test-beam conditions, the CAL single-particle relative energy resolutions were 
a{E)/E = 0.18/v 7 ^ for electrons and cr{E)/E = 0.35/v 7 ^ for hadrons, with E in GeV. 

The luminosity was measured from the rate of the bremsstrahlung process ep — > e^p. The 
resulting small-angle energetic photons were measured by the luminosity monitor [16], a 
lead-scintillator calorimeter placed in the HERA tunnel at Z = — 107 m. 
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4 Data selection 



The data were collected during the running period 1998-2000, when HERA operated with 
protons of energy E p = 920 GeV and electrons or positrons 2 of energy E e = 27.5 GeV, 
and correspond to an integrated luminosity of 81.7 ± 1.9 pb _1 . 

Neutral current DIS events were selected offline using criteria similar to those reported 
previously [5]. The main steps are given below. 

A reconstructed event vertex consistent with the nominal interaction position was re- 
quired and cuts based on tracking information were applied to reduce the contamination 
from beam-induced and cosmic-ray background. The scattered-electron candidate was 
identified using the pattern of energy deposits in the CAL [17]. The energy, E' e , and polar 
angle, 6 e , of the electron candidate were also determined from the CAL measurements. 
The double-angle method [18], which uses 6 e and an angle 7 that corresponds, in the 
quark-parton model, to the direction of the scattered quark, was used to reconstruct Q 2 . 
The angle 7 was reconstructed using the CAL measurements of the hadronic final state. 

Electron candidates were required to have an energy E' e > 10 GeV, to ensure a high and 
well understood electron-finding efficiency and to suppress background from photopro- 
duction. The inelasticity variable, y, as reconstructed using the electron energy and polar 
angle, was required to be below 0.95; this condition removed events in which fake electron 
candidates from photoproduction background were found in the FCAL. The requirement 
38 < (E — pz) < 65 GeV, where E is the total CAL energy and pz is the Z compo- 
nent of the energy measured in the CAL cells, was applied to remove events with large 
initial-state radiation and to reduce further the photoproduction background. Remaining 
cosmic rays and beam-related background were rejected by requiring the total missing 
transverse momentum, p™ lss , to be small compared to the total transverse energy, E^, 
PT i88 /\/Er t < 3 VGeV. The kinematic range was restricted to Q 2 > 125 GeV 2 . 

The kx cluster algorithm was used in the longitudinally invariant inclusive mode to recon- 
struct jets in the measured hadronic final state from the energy deposits in the CAL cells. 
The jet algorithm was applied after excluding those cells associated with the scattered- 
electron candidate. Jet transverse-energy corrections were computed using the method 
developed in a previous analysis [5]. Events were required to have at least one jet of 
E^ > 14 GeV and -1 < r/ jet < 2.5. The final sample of 128986 events contained 132818 
jets, of which 21162 jets had exactly two subjets at y cut = 0.05. 



2 In the following, the term "electron" denotes generically both the electron (e ) and the positron (e + ). 
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5 Monte Carlo simulation 



Samples of events were generated to determine the response of the detector to jets of 
hadrons and the correction factors necessary to obtain the hadron-level subjet cross 
sections. The hadron level is defined as those hadrons with lifetime r > 10 ps. The 
generated events were passed through the Geant 3.13-based [19] ZEUS detector- and 
trigger-simulation programs [13]. They were reconstructed and analysed applying the 
same program chain as to the data. 

Neutral current DIS events including radiative effects were simulated using the Her- 
acles 4.6.1 [20] program with the Djangoh 1.1 [21] interface to the hadronisation 
programs. Heracles includes corrections for initial- and final-state radiation, vertex 
and propagator terms, and two-boson exchange. The QCD cascade is simulated using 
the colour-dipole model (CDM) [22] including the LO QCD diagrams as implemented in 
Ariadne 4.08 [23] and, alternatively, with the MEPS model of Lepto 6.5 [24]. The 
CTEQ5D [25] proton PDFs were used for these simulations. Fragmentation into hadrons 
is performed using the Lund string model [26] as implemented in Jetset [27,28]. 

The jet search was performed on the Monte Carlo (MC) events using the energy measured 
in the CAL cells in the same way as for the data. The same jet algorithm was also applied 
to the final-state particles (hadron level) and to the partons available after the parton 
shower (parton level) to compute hadronisation correction factors (see Section 6). 

6 QCD calculations 

The 0(a 2 s ) NLO QCD calculations used to compare with the data are based on the 
program Disent [29]. The calculations used a generalised version of the subtraction 
method [30] and were performed in the massless MS renormalisation and factorisation 
schemes. The number of flavours was set to five; the renormalisation and factorisation 
(Hf) scales were set to fiR — fiF — Q] ol s was calculated at two loops using A|| = 220 MeV 
which corresponds to a s (Mz) = 0.118. The ZEUS-S [31] parameterisations of the proton 
PDFs were used. The results obtained with Disent were cross-checked by using the 
program NLOJET++ [32]. 

Since the measurements refer to jets of hadrons, whereas the QCD calculations refer to 
jets of partons, the predictions were corrected to the hadron level using the MC sam- 
ples described in Section 5. The multiplicative correction factor, Chad, defined as the 
ratio of the cross section for subjets of hadrons to that of partons, was estimated with 
the Lepto-MEPS model, since it reproduced the shape of the QCD calculations better. 
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The normalised cross-section calculations changed typically by less than ±20% upon ap- 
plication of the parton-to-hadron corrections, except at the edges of the distributions, 
where they changed by up to ±50%. Other effects not accounted for in the calculations, 
namely QED radiative corrections and Z° exchange, were found to be very small for the 
normalised cross-section calculations and neglected. 

The following theoretical uncertainties were considered (as examples of the size of the 
uncertainties, average values of the effect of each uncertainty on the normalised cross 
section as functions of E^/E^, r] sh i —rf ct , \<p sb ^ — <ft ct \ and a sh i are given in parentheses): 

• the uncertainty in the modelling of the parton shower was estimated by using different 
models (see Section 5) to calculate the parton-to-hadron correction factors (5.6%, 
13.2%, 7.6%, 5.3%); 

• the uncertainty on the calculations due to higher-order terms was estimated by varying 
fi R by a factor of two up and down (0.01%, 0.46%, 0.58%, 0.34%); 

• the uncertainty on the calculations due to the choice of [ip was estimated by varying 
Hp by a factor of two up and down (0.05%, 0.43%, 0.11%, 0.12%); 

• the uncertainty on the calculations due to those on the proton PDFs was estimated 
by repeating the calculations using 22 additional sets from the ZEUS analysis [31]; 
this analysis takes into account the statistical and correlated systematic experimental 
uncertainties of each data set used in the determination of the proton PDFs (0.07%, 
0.18%, 0.12%, 0.05%); 

• the uncertainty on the calculations due to that on a s [Mz) was estimated by repeating 
the calculations using two additional sets of proton PDFs, for which different values 
of a s (Mz) were assumed in the fits. The difference between the calculations using 
these various sets was scaled to reflect the uncertainty on the current world average 
of a, [33] (0.02%, 0.04%, 0.05%, 0.01%). 

These uncertainties were added in quadrature and are shown as hatched bands in the 
figures. 

7 Corrections and systematic uncertainties 

The sample of events generated with CDM, after applying the same offline selection as 
for the data, gives a reasonably good description of the measured distributions of the 
kinematic, jet and subjet variables; the description provided by the MEPS sample is 
somewhat poorer. The comparison of the measured subjet distributions and the MC 
simulations is shown in Fig. 2. 
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The normalised differential cross sections were obtained from the data using the bin-by-bin 
correction method, 

adA~ aC-AAi' ' 

where iVdata,; is the number of subjets in data in bin % of the subjet variable A, N^q a (N^c i) 
is the number of subjets in MC at hadron (detector) level, C is the integrated luminosity 
and AAi is the bin width. The MC samples of CDM and MEPS were used to compute 
the acceptance correction factors to the subjet distributions. These correction factors 
took into account the efficiency of the trigger, the selection criteria and the purity and 
efficiency of the jet and subjet reconstruction. 

The following sources of systematic uncertainty were considered for the measured subjet 
cross sections (as examples of the size of the uncertainties, average values of the effect 
of each uncertainty on the normalised cross section as functions of E^ s /E^, r] sh i — rf et , 
|0sbj _ 0jct| an( j ^sbj are gj ven j n parentheses): 

• the deviations in the results obtained by using either CDM or MEPS to correct the 
data from their average were taken to represent systematic uncertainties due to the 
modelling of the parton shower (0.5%, 2.9%, 2.6%, 1.3%); 

• variations in the simulation of the CAL response to low-energy particles (0.3%, 1.6%, 
1.2%, 0.6%). 

Other uncertainties, such as those arising from the uncertainty in the absolute energy scale 
of the jets [1,34], the uncertainty in the simulation of the trigger and the uncertainty in 
the absolute energy scale of the electron candidate [35], were investigated and found to 
be negligible. The systematic uncertainties were added in quadrature to the statistical 
uncertainties and are shown as error bars in the figures. 



8 Results 

Normalised differential subjet cross sections were measured for Q 2 > 125 GeV 2 for jets 
with > 14 GeV and -1 < rf ct < 2.5 which have exactly two subjets for y cnt = 0.05. 

The distribution of the fraction of transverse energy, (1 / 'a) (da j 'd(E^ / 'E^)) , is presented 
in Fig. 3a. It contains two entries per jet and is symmetric with respect to E^/E^ = 0.5 
by construction. This distribution has a peak for 0.4 < E^/E^ < 0.6, which shows that 
the two subjets tend to have similar transverse energies. 
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The ?7 sb j — t]^ ct data distribution is shown in Fig. 3b and also has two entries per jet. The 
measured cross section has a two-peak structure; the dip around r] sb ^ — rf ct = is due to 
the fact that the two subjets are not resolved when they are too close together. 

Figure 3c presents the measured normalised cross section as a function of |0 sb j — ^ et |. 
There are two entries per jet in this distribution. The distribution has a peak for 0.2 < 
|0sbj _ 0jet| < q 3. foe suppression around \<j) sh ^ — <ft et \ = also arises from the fact that 
the two subjets are not resolved when they are too close together. 

The data distribution as a function of a sh i (one entry per jet) increases as a * increases 
(see Fig. 3d). This shows that the subjet with higher transverse energy tends to be 
in the rear direction. This is consistent with the asymmetric peaks observed in the 
^sbj _ ^jct distribution (see Fig. 3b). Figure 4 shows the 77 j — rf ct distribution for those 
jets which have two subjets with asymmetric E^ (E^i ow / E^ 1 < 0.4, or, equivalently, 
-^high/^T* > 0-6)) separately for the subjet with higher and lower E^ 1 . It is to be noted 
that since the jet axis is reconstructed as the transverse-energy-weighted average of the 
subjet axes, the subjet with higher E^ is constrained to be closer to the jet axis than 
that of the lower E^ subjet. The measured distributions show that the higher (lower) 
E^ subjet tends to be in the rear (forward) direction. All these observations support the 
expectation of the presence of colour- coherence effects between the initial and final states 
and, in particular, the tendency of the subjet with lower E^ 1 to be emitted predominantly 
towards the proton beam direction. 

8.1 Comparison with NLO QCD calculations 

Next-to-leading-order QCD calculations are compared to the data in Figs. 3 and 4. The 
QCD predictions give an adequate description of the data. However, the data points are 
situated at the upper (lower) edge of the theoretical uncertainty in some regions of the 
subjet variables such as E^^/E^ ~ 0.5, \<fi sh i — <^ ct | ~ 0, a sh i ~ and the peaks in 
the r/ sbj - rf ct distribution {EpjEjf ~ 0.25, |0 sbj - ft ct \ > 0.3 and |r/ sb J - rf ct \ > 0.5). 
Since the calculations are normalised to unity, the uncertainties are correlated among the 
points; this correlation gives rise to the pulsating pattern exhibited by the theoretical 
uncertainties. 

The calculation of the cross section as a function of E^/E^ 1 exhibits a peak at 0.4 < 
E^ j Eq? < 0.6, as seen in the data. The calculations for the r] sb ^ — rf 01 and a sh i dis- 
tributions predict that the subjet with higher transverse energy tends to be in the rear 
direction, in agreement with the data. This shows that the mechanism driving the subjet 
topology in the data is the eq — > eqg and eg —* eqq subprocesses as implemented in the 
pQCD calculations. 
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To gain further insight into the pattern of parton radiation, the predictions for quark- and 
gluon-induced processes (see Section 2) are compared separately with the data in Fig. 5. 
The NLO calculations predict that the two-subjet rate is dominated by quark-induced 
processes: the relative contribution of quark- (gluon-) induced processes is 81% (19%). 
The shape of the predictions for these two types of processes are different; in quark-induced 
processes, the two subjets have more similar transverse energies (see Fig. 5a) and are closer 
to each other (see Fig. 5b and 5c) than in gluon-induced processes. The comparison with 
the measurements shows that the data are better described by the calculations for jets 
arising from a qg pair than those coming from a qq pair. 



8.2 , Q 2 and x dependence of the subjet distributions 

Figures 6 to 9 show the normalised differential subjet cross sections in different regions 
of Ej,. Even though the mean subjet multiplicity decreases with increasing E^ [4], the 
measured normalised differential subjet cross sections have very similar shapes in all E^ 
regions for all the observables considered. This means that the subjet topology does not 
change significantly with Ej, . This is better illustrated in Fig. 10, where the data for all 
Eji regions are plotted together. In particular, it is observed that the maximum of each 
measured normalised cross section in every region of E^ occurs in the same bin of the 
distribution. To quantify the Ej, dependence more precisely, Fig. 11 shows the maximum 
value of the measured normalised cross section for each observable as a function of E^ 
together with the NLO predictions. The spread of the measured maximum values of the 
normalised cross sections is ±(4 — 6)%. For each observable, the scaling behaviour of the 
normalised differential subjet cross sections is clearly observed and in agreement with the 
expectation that the splitting functions depend weakly on the energy scale. The NLO 
QCD calculations are in agreement with the data and support this observation. 

Figures 12 to 15 show the normalised differential subjet cross sections in different regions 
of Q 2 . In this case, it is observed that while the shape of the E^/E^ distribution does not 
change significantly with Q 2 , some dependence can be seen in the other observables. For 
example, the dip in the r] sh] - — if et distribution is shallower for 125 < Q 2 < 250 GeV 2 than 
at higher Q 2 and the shape of the a sh] distribution for 125 < Q 2 < 250 GeV 2 is somewhat 
different than for the other regions (see Fig. 16). These features of the data are reasonably 
reproduced by the NLO QCD calculations and understood as a combination of two effects: 
the fraction of gluon-induced events is predicted to be 32% for 125 < Q 2 < 250 GeV 2 
and below 14% for higher Q 2 ; the shape of the normalised cross sections as functions of 
^sbj _ r/j et and a sbj changes from the region 125 < Q 2 < 250 GeV 2 to 250 < Q 2 < 500 GeV 2 

(see Fig. 17) for quark- and gluon-induced events. It is observed that the maximum of 
each measured normalised cross section in every region of Q 2 occurs in the same bin 
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of the distribution, except for \<p sh i — <ft ct \ in the highest-Q 2 region. Figure 18 shows the 
maximum 3 value of the measured normalised cross section for each observable as a function 
of Q 2 together with the NLO predictions. The spread of the measured maximum values 
of the normalised cross sections as functions of Ep / and |0 sbj - <^ ct | is ±(3 - 4)%. 
On the other hand, the measured and predicted maximum values for the normalised cross 
sections as functions of rf * — rf et and a sbl exhibit a step-like behaviour between the 
lowest-Q 2 region and the rest. 

Figures 19 to 22 show the normalised differential subjet cross sections in different regions 
of x. Figure 23 shows the data for all x regions plotted together. It is observed that the 
maximum of each measured normalised cross section in every region of x occurs in the 
same bin of the distribution, except for \<p sh> — <ft ct \ in the highest x region. Figure 24 
shows the maximum 3 value of the measured normalised cross section for each observable 
as a function of x. The shape of the Ep/E^ measured distribution does not change 
significantly with x, whereas some dependence is expected (see Fig. 24a). The dependence 
of the ?7 sb j — rf ct and a sh i distributions with x exhibits features similar to those observed 
in the study of the Q 2 dependence; in particular, the maximum values (see Figs. 24b and 
24d) exhibit a monotonic increase as x increases, which is reasonably reproduced by the 
calculations. As discussed previously, these features are understood as a combination of 
two effects: a decrease of the predicted fraction of gluon-induced events from 44% for 
0.004 < x < 0.009 to 6% for x > 0.093 and the change in shape of the normalised cross 
sections for quark- and gluon-induced processes as x increases (see Fig. 25). 

To investigate the origin of the change in shape of the normalised differential cross sections 
between the lowest and higher Q 2 and x regions, LO and NLO calculations were compared. 
The most dramatic change is observed when restricting the kinematic region to 125 < 
Q 2 < 250 GeV 2 or 0.004 < x < 0.009 (see Fig. 26); the LO calculation of the r/ sbj - r/ jct 
distribution does not exhibit a two-peak structure as seen in the NLO prediction and in 
the data. In addition, the LO calculation of the a sh i distribution peaks at a sh i ~ ir/2 in 
contrast with the NLO prediction and the data. This proves that the NLO QCD radiative 
corrections are responsible for these variations in shape and necessary for describing the 
data. 

In summary, while the shapes of the normalised differential cross sections show only a weak 
dependence on Ej. , their dependence on Q 2 and x have some prominent features at low Q 2 
or x. The weak dependence on Ej? is consistent with the expected scaling behaviour of the 
splitting functions; however, the restriction to low Q 2 or x values demonstrates that the 
NLO QCD radiative corrections are important there. The NLO QCD calculations, which 
include the two competing processes eq — > eqg and eg — > eqq and radiative corrections, 

3 For the \cf> sh ' 3 — <ft et \ distribution, the same bin has been used for consistency. 
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adequately reproduce the measurements. 



9 Summary 

Normalised differential subjet cross sections in inclusive-jet NC DIS were measured in 
ep collisions using 81.7 pb _1 of data collected with the ZEUS detector at HERA. The 
cross sections refer to jets identified in the laboratory frame with the kr cluster algorithm 
in the longitudinally invariant inclusive mode and selected with > 14 GeV and 
— 1 < rf ct < 2.5. The measurements were made for those jets which have exactly two 
subjets for y cnt = 0.05 in the kinematic region defined by Q 2 > 125 GeV 2 . 

The cross sections were measured as functions of E^/E^, r] sh i — if et , \<jf h ' s — <fi iCt \ and a shi . 
The data show that the two subjets tend to have similar transverse energies and that the 
subjet with higher transverse energy tends to be in the rear direction. This is consistent 
with the effects of colour coherence between the initial and final states, which predict that 
soft parton radiation is emitted predominantly towards the proton beam direction. 

An adequate description of the data is given by NLO QCD calculations. This means that 
the pattern of parton radiation as predicted by QCD reproduces the subjet topology in 
the data. Furthermore, the subjet distributions in the data are better described by the 
calculations for jets arising from a quark-gluon pair. 

The normalised cross sections show a weak dependence on E^ , in agreement with the 
expected scaling behaviour of the splitting functions. By restricting the measurements to 
low Q 2 or x values, significant differences in shape are observed, which can be primarily 
attributed to NLO QCD radiative corrections. 
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Figure 2: Detector-level normalised subjet data distributions (dots) for jets with 
&f > U GeV and -1 < t/J * < 2.5 which have two subjets for y cnt = 0.05 in 
the kinematic region given by Q 2 > 125 GeV 2 as functions of (a) Ep / , (b) 
rf° 1 —rf et , (c) |0 sb j — 0j ct | and (d) a sh K The statistical uncertainties are smaller than 
the marker size. For comparison, the distributions of the CDM (solid histograms) 
and MEPS (dot-dashed histograms) Monte Carlo models are included. 
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Figure 3: Measured normalised differential subjet cross sections (dots) for jets 
with > 14 GeV and —1 < rf ct < 2.5 which have two subjets for y cut = 0.05 
in the kinematic region given by Q 2 > 125 GeV 2 as functions of (a) E^ s / Erf? , 
(b) r] sh i — rf ct , (c) |0 1 — <^ ct | and (d) a sh K The inner error bars represent the 
statistical uncertainties of the data, the outer error bars show the statistical and 
systematic uncertainties added in quadrature. In many cases, the error bars are 
smaller than the marker size and are therefore not visible. For comparison, the 
NLO QCD predictions (solid histograms) are included. The hatched bands represent 
the theoretical uncertainty. 
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Figure 5: Measured normalised differential subjet cross sections (dots) for jets 

with > 14 GeV and —1 < rf ct < 2.5 which have two subjets for y cut = 0.05 in 

the kinematic region given by Q 2 > 125 GeV 2 as functions of (a) Ep/E^, (b) 
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Figure 6: Measured normalised differential subjet cross sections (dots) for jets 
with E^ > 14 GeV and -1 < rf ct < 2.5 which have two subjets for y cnt = 0.05 in 
the kinematic region given by Q 2 > 125 GeV 2 as functions of E^ / E^ in different 
regions of E^ . Other details are as in the caption to Fig. 3. 
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with £p > 14 GeFand -1 < 77-)°* < 2.5 which have two subjets for y CVLt = 0.05 in 
the kinematic region given by Q 2 > 125 GeV 2 as functions of r] sh i — rf ct in different 
regions of ■ Other details are as in the caption to Fig. 3. 
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Figure 8: Measured normalised differential subjet cross sections (dots) for jets 
with E!£ > 14 GeV and — 1 < rf ct < 2.5 which have two subjets for y cut = 0.05 
in the kinematic region given by Q 2 > 125 GeV 2 as functions of \<f) sh i — <ft et \ in 
different regions of ■ Other details are as in the caption to Fig. 3. 
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Figure 9: Measured normalised differential subjet cross sections (dots) for jets 
with > 14 GeV and -1 < T]i ct < 2.5 which have two subjets for y cut = 0.05 
in the kinematic region given by Q 2 > 125 GeV 2 as functions of a sh> in different 
regions of . Other details are as in the caption to Fig. 3. 
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Figure 10: Measured normalised differential subjet cross sections for jets with 
> 14 GeV and — 1 < rf ct < 2.5 which have two subjets for y cut = 0.05 in 
the kinematic region given by Q 2 > 125 GeV 2 as functions of (a) / , (b) 
rj sh ^—rj^ ct , (c) \(jf hi — 0j ct | and (d) a shs in different regions of E^ . Details concerning 
the error bars are as in the caption to Fig. 3. 
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Figure 11: Maximum of the measured normalised differential (a) E^ 3 / ' Ej, , (b) 
^sbj _ ^jct^ |0sbj _ 0jet| an( ^ a sb] su f)j e f cross sections (dots) for jets with 

E^ > 14 GeV and -1 < rf ct < 2.5 which have two subjets for y cut = 0.05 in the 
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histograms) processes are also shown separately. Other details are as in the caption 
to Fig. 3. 
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Figure 12: Measured normalised differential subjet cross sections (dots) for jets 
with > 14 GeV and — 1 < if et < 2.5 which have two subjets for y cut = 0.05 in 
the kinematic region given by Q 2 > 125 GeV 2 as functions ofE^/E^ in different 
regions of Q 2 . Other details are as in the caption to Fig. 3. 
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Figure 13: Measured normalised differential subjet cross sections (dots) for jets 
with E 3 ^ > 14 GeV and — 1 < rj :,ct < 2.5 which have two subjets for y cut = 0.05 in 
the kinematic region given by Q 2 > 125 GeV 2 as functions of r] sh] - — rf Qt in different 
regions of Q 2 . Other details are as in the caption to Fig. 3. 
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Figure 14: Measured normalised differential subjet cross sections (dots) for jets 
with > 14 GeV and —1 < rf ct < 2.5 which have two subjets for y cut = 0.05 
in the kinematic region given by Q 2 > 125 GeV 2 as functions of \<j) sh i — <ft ct \ in 
different regions of Q 2 . Other details are as in the caption to Fig. 3. 
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Figure 15: Measured normalised differential subjet cross sections (dots) for jets 
with E 3 rf l > 14 GeV and — 1 < if et < 2.5 which have two subjets for y cut = 0.05 
in the kinematic region given by Q 2 > 125 GeV 2 as functions of a sh> in different 
regions of Q 2 . Other details are as in the caption to Fig. 3. 
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Figure 16: Measured normalised differential subjet cross sections for jets with 
EP* > 14 GeV and -1 < rf ct < 2.5 which have two subjets for y cut = 0.05 in 

the kinematic region given by Q 2 > 125 GeV 2 as functions of (a) Ep/E^\ (b) 
^sbj_^jet^ |0sbj_ ( ^jet| an ^ ^sbj ^ n di^ eren i regions of Q 2 . Details concerning 

the error bars are as in the caption to Fig. 3. 
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Figure 17: Predicted normalised differential subjet cross sections (solid his- 
tograms) for jets with > 14 GeV and —1 < r]^ ct < 2.5 which have two subjets 
for y cut = 0.05 in the kinematic region given by Q 2 > 125 GeV 2 as functions of 
(a,c) r] sb i — rf ot and (b,d) a sh] in different regions of Q 2 . The NLO predictions for 
quark- (dotted histograms) and gluon-induced (dot-dashed histograms) processes are 
also shown separately. 
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Figure 18: Maximum of the measured normalised differential (a) E^ i / E^ , (b) 
r) i — rf ct , (c) \(jf h ' s — (ft ct \ and (d) a shs subjet cross sections (dots) for jets with 
E^ > 14 GeV and -1 < rf ct < 2.5 which have two subjets for y cut = 0.05 in the 
kinematic region given by Q 2 > 125 GeV 2 as a function of Q 2 . For comparison, 
the NLO predictions for quark- (dotted histograms) and gluon-induced (dot-dashed 
histograms) processes are also shown separately. Other details are as in the caption 
to Fig. 3. 
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Figure 19: Measured normalised differential subjet cross sections (dots) for jets 
with > 14 GeV and — 1 < rf et < 2.5 which have two subjets for y cut = 0.05 in 
the kinematic region given by Q 2 > 125 GeV 2 as functions of E^/E^ in different 
regions of x. Other details are as in the caption to Fig. 3. 
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Figure 20: Measured normalised differential subjet cross sections (dots) for jets 
with > 14 GeV and — 1 < if et < 2.5 which have two subjets for y cut = 0.05 in 
the kinematic region given by Q 2 > 125 GeV 2 as functions of r] sh i — rf ct in different 
regions of x. Other details are as in the caption to Fig. 3. 
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Figure 21: Measured normalised differential subjet cross sections (dots) for jets 
with E 3 rf l > 14 GeV and — 1 < rf et < 2.5 which have two subjets for y cut = 0.05 
in the kinematic region given by Q 2 > 125 GeV 2 as functions of \<f) sh i — <ft et \ in 
different regions of x. Other details are as in the caption to Fig. 3. 



35 



ZEUS 



|0.5 

D 
T3 




0.5 



0.004 < x < 0.009 



/ ml I I I 



mi 



mm 



0.009 < x < 0.019 



0.019 < x < 0.043 



0.043 < x < 0.093 



I I I LI I 17 



1 1 1 a z z z 



W7^ 




0.5 



0.093 < x < 0.447 

• ZEUS 82 pb 1 

WFFFR NLO 



"7f7~' 







oc sbj (rad) 



a sbj (rad) 



Figure 22: Measured normalised differential subjet cross sections (dots) for jets 
with > 14 GeV and — 1 < rf ct < 2.5 which have two subjets for y cut = 0.05 
in the kinematic region given by Q 2 > 125 GeV 2 as functions of a sh i in different 
regions of x. Other details are as in the caption to Fig. 3. 
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Figure 23: Measured normalised differential subjet cross sections for jets with 
E^ > 14 GeV and -1 < rf ct < 2.5 which have two subjets for y cut = 0.05 in 

the kinematic region given by Q 2 > 125 GeV 2 as functions of (a) E^ s / , (b) 
^sbj _^jct^ |0sbj _0jct| an ^ ^ ^sbj ^ n different regions of x. Details concerning 
the error bars are as in the caption to Fig. 3. 
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Figure 24: Maximum of the measured normalised differential (a) E^ i / E^ , (b) 
r) i — rf ct , (c) |0 j — (^ et | and ^ a sb -j subjet cross sections (dots) for jets with 
E^ > 14 GeFand -1 < rf et < 2.5 which have two subjets for y cut = 0.05 in the 
kinematic region given by Q 2 > 125 GeV 2 as a function of x. For comparison, 
the NLO predictions for quark- (dotted histograms) and gluon-induced (dot-dashed 
histograms) processes are also shown separately. Other details are as in the caption 
to Fig. 3. 
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Figure 25: Predicted normalised differential subjet cross sections (solid his- 
tograms) for jets with > 14 GeV and — 1 < if et < 2.5 which have two subjets 
for y CVLt = 0.05 in the kinematic region given by Q 2 > 125 GeV 2 as functions of 
(a,c) r] ] — rf ct and (b,d) a sb i in different regions of x. The NLO predictions for 
quark- (dotted histograms) and gluon-induced (dot-dashed histograms) processes are 
also shown separately. 
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Figure 26: Measured normalised differential subjet cross sections (dots) for jets 
with > 14 GeV and — 1 < rf et < 2.5 which have two subjets for y cut = 0.05 
in restricted Q 2 and x regions as functions of (a,c) 7] shs — rf ct and (b,d) a sb K The 
NLO (solid histograms) and LO (dashed histograms) calculations are also shown. 
The hatched bands represent the NLO theoretical uncertainty. 
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